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Abstract 

Call a spectrum of Hamiltonian sparse if each eigenvalue can be quickly restored with 
accuracy e from its rough approximation in within ei by means of some classical algorithm. 
It is shown how a behavior of system with sparse spectrum up to time T = can be 
predicted with fidelity p on quantum computer in time t = plus the time of 

classical algorithm. The quantum knowledge of Hamiltonian H eigenvalues is considered 
as a wizard Hamiltonian Wh which action on any eigenvector of H gives the corresponding 
eigenvalue. Speedup of evolution for systems with sparse spectrum is possible because for 
such systems wizard Hamiltonians can be quickly simulated on a quantum computer. This 
simulation, generalizing Shor trick, is a part of presented algorithm. In general case the 
action of wizard Hamiltonian cannot be simulated in time smaller than the dimension of 
main space which is exponential of the size of quantum system. For an arbitrary system 
(even for classical) its behavior cannot be predicted on quantum computer even for one 
step ahead. This method can be used also for restoration of a state of an arbitrary primary 
system in time instant —T in the past with the same fidelity which requires the same time. 



1 Introduction and background 

A behavior of typical quantum system cannot be analyzed at hand or by classical computer be- 
cause huge dimension of its Hamiltonian makes it impossible to solve Shroedinger equation even 
numerically. Nevertheless one can expect that this analysis would be easier in the framework 
of quantum computing. The idea of this approach is to force one quantum system to simulate 
a behavior of other more simple primary system with some profit in time. In the particular 
case where the primary system is a classical computer with oracle this is a problem of quantum 
speedup for classical computations: given a classical algorithm with oracle, is there a quantum 
algorithm computing the same function faster using the same oracle quantumly? Examples: 
Shor factoring algorithm (look at Grover search algorithm (look at |(?r|). 



In general case the problem of prediction has the following form: given a low of evolution 
for some system (classical or quantum) is there a device (wizard) of not exponential size which 
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can predict the behavior of this system? To predict the behavior means that given a time 
instant, wizard returns a state of initial system in this instant earher than this state appears 
in its natural evolution. Note that in general case of quantum evolution a wizard may use 
quantumly not only Hamiltonian (like in case of quantum computations with oracle) but some 
hidden information about primary system. What kind of such information may be useful for 
prediction? This is the information about eigenvalues of primary Hamiltonian. This article 
analyses the possibility of quantum predictions in terms of spectrum of primary system. 

Quantum method of finding eigenvalues was first presented by Abrams and Lloyd in their 
work ||AL|| . Their method, generalizing a quantum part of Shor factoring algorithm, uses Quan- 
tum Fourier transform (QFT) and requires of order N implementations of initial Hamiltonian 
where is the dimension of the main space. Thus, their algorithm requires exponential time 
and it cannot be used for quantum speedup. However, for the systems whose spectrum is 
sparse the idea of this approach with QFT may be used for predictions. I proceed with the 
exact definitions. 



Let H be some Hamiltonian which induces unitary transformations U in some A^ dimensional 
space Y of states of n qubits by the usual rule in quantum mechanics: U{t) = exp^—iHr/h), 
N = 2"'. In order to apply notions from the algorithm theory, such as complexity, etc., we 
assume that the time of Hamiltonian action is always r = 1. We may of course choose another 
value for r but this would only involve the change in time scale without any subsequences 
for the algorithm complexity. Call a system primary if its evolution is determined by unitary 
transformation U. 

We assume that the complexity of computation is the number of U applications. Time of all 
other transformations we use here is assumed to be negligible comparatively with the time of 
computation. Quantum Fourier transform (QFT) requires the time of order (look at |[Sh|| ). 
Here we shall use this transform for log2(M) gubits in the following form: 



Let eigenvalues of H have the form —2'KUJkh , /c = 0, 1, . . . , A^ — 1, when frequencies ujk are 
real numbers from the segment [0, 1). Then the eigenvalues of U will be exp(27ricc;fc). This is 
not loss in generality because we always can choose other unit r for Hamiltonian action time. 
Let us give the precise definitions of the "knowledge of eigenvalues" of H . It means that there 
exists another Hamiltonian Wh called a wizard for H which acts in A^^ dimensional Hilbert 
space and returns an approximation of frequencies uJk in within 1/A^ given eigenvector $fc of 
H . Introduce the following notations. Every frequency uJk has binary notation 0.eie2 . . .. Let 
p < nhe some integer, M = T' . Denote the string of ones and zeroes t\e2 . . . Cp by . may 
be considered as integer if we cancel all first zeroes. And vise versa, every integer / less than 
M can be written in form if we add suitable number of zeroes in front. Then the number 



2 Definitions 
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uj^ = 0.e| will be approximation of Wk in within 1/M. We shall define this number O.e^ from 
[0, 1) by {0.l)p. Wizard action in these notations will be 

Wu\^k,b)^\^k,b^el) 

where is the biwise addition modulo 2. 

At last introduce auxiliary transformations WH and U^g^. Let the memory is divided into 
the main part a; of n qubits and ancilla a of p < n qubits: \x,a), M = 2^. Put: 

M-l 

1) . WH|x, a) = -4^ J2 (— l)" ''|a;, s). This is Walsh - Hadamard transform applied to ancilla 

where • denotes dot product modulo 2. 

2) . l]P^g\x,a) = |U"a;, a). This is the result of a sequential applications of U to the main 
register. 

Walsh-Hadamard transform can be fulfilled in a standard model of quantum computer. 
To fulfill Ufgg it would suffice to use the following oracle Ucond (conditional application of U) 
depending on U: Ucond\x,a) — > l^^',^), where x' = x if a = and x' = Ux if o; = 1, a is 
one qubit register. An application of Ucond cannot be reduced to the simple using of U as an 
oracle because a conditional application of U is quantumly controlled by the second register. 
One proposal about its practical implementation can be found in the section 4.5. 



3 Wizard transformation 



3.1 How wizard predicts evolution 

Assume temporarily that we have exact equations cuk — oJk- Given a wizard transformation 
Wu how can we predict evolution of initial system? Let |^) denote initial state as a contents of 
n qubits main register. Let 

N-1 

^^Y.^k\^k (1) 

A:=0 

be the expansion of our state in basis of eigenvectors of U . Add n qubits ancillary register 
initiafized by zeroes and obtain the state 0"). 

Now apply the wizard transformation Wu to the main register {p — n). It gives the state 

Af-l 
A;=0 

Given the number t we can turn each state in the first register by angle determined by the 
second register and t: 27ru;kt. This gives the state 

N-1 

i" = J2 Xkexp{2niujkt)\^k,el). 

k=0 

At last apply wizard again obtaining 

N-1 

C= E^A^exp(27riu;fei)|$fc,0"), 

fc=0 
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and now cancellation ancilla gives exactly the state U^\C,) which is the state in time instant t 
of initial system. Thus given a wizard we can predict the behavior of initial system provided 
the action of wizard takes smaller time than t/2. 



3.2 Simulation a wizard 

In this subsection I describe the procedure of simulation a wizard action on zero ancilla. The 
particular case p = n of this procedure was proposed in other purpose in the paper |[AL| . 



We start from the state of the form (1) with ancilla attached initialized by zeroes. Wizard 
action on the initial state with zero ancilla is defined by the following: 

QFT,, UL,WH|e,0^). (2) 

Why this procedure must work? Let us again assume temporarily that Uk = Cj^^ ex- 

N-lM-l 

actly. Then we have WH|^,0^) = Xo = -jti Xk\^k-,s). Application of Ufg^ gives: 



fc=0 s=0 

Af-lM-l Af-1 

Xi = UfecXo = 7^ E E a;fcexp(27riiu;fcs)|$fc, s). Now QFTj^^xi = E Xk\^k,el). This is just 

the wizard action. But in general case the result of QFT transform has not so simple form. It 
means that we cannot merely clean ancilla by repetition of wizard action in predicting algo- 
rithm even if we have a wizard action on all words. Hence the precision of a wizard simulation 
by (2) must be elaborated in more details. 

3.3 Accuracy of a wizard simulation 

Let {Cjk,i} be some set of integers. Denote Le{u)k,i) = {i '■ \{0.ijk,i)p — ujk\ < e or 
\{0.ujk,i)p — uJk — l\ < £}■ Let be a state of the form (1). 
Definition A transformation W of the form 

N-l M-1 

k=0 i=0 

N-l 

is called a transformation of Ws^e type if J2 J2 l-^i.fcP > 1 — S for any ^. Thus, 

6 is an error probability from the quantum superposition, and £ is a precision of eigenvalues 
approximations. 

Lemma 1 The transformation QFT^ U^g^ WH belongs to the type W± k_ for any K e 
{1,2,. ..,M}. 

Proof 

Denote QFT^ U^^^ WH |^, 0^) by Xp,2- Then we have 

N-lM-l 

Xp,2 = V7 E E XkHi^k\^k,l) (3) 
k=0 1=0 
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M-1 

where Hik— J2 exp(27ris(a;fe — (OJ)p)). Put A = a;^ — (OJ)p. Then by summing the progression 

s=0 

we obtain 

_ 1 - exp(27r^MA) 
- l-cxp(27r^A) ■ 

Let A = ^ + 5, K = K{1), where < 5 < 1/M, K integer. Then each Hi ^ depends on 
K : k = Hi^k{K). Fix some value for K: Kq. It means that the accuracy of eigenvalues 
approximations will be Kq/M. Estimate the sum of squared amplitudes for states |$fc,0 ^^r 
which K : Kq < K < M — Kq. Using that module of denominator in (4) is separated from 
zero by min{^, ^^^^}, we have: 

iV-l M-Ko 

E E M^\Hi,k{K)\'' < 

k=0 K=Ko 

N-IM-Ka M-Ko 

4 E E IxkWn'K' < ^ E w<h 

k=0 K=Ko K=Ko ° 

This is exactly what is needed. Lemma 1 is proved. 

What would happen if we apply other sequence of transformations QFT^^ (^^eq)~^ 
instead of QFTji^ U^^^ WH for reveahng frequencies? 

Lemma 2 Let for some state x 

M-l 

X2 = QFT^ U?,, WH|x, (F) = E E VkA^k, I) 

k 1=0 

M-l 

=QFT^i (u^J-i WH|x,o^') =E E yii\^k,l) 



k 1=0 



PutSk,i^uj^-{0.l)p,A^u;k-{0.l)p,X2 = J: E Z/;iexp(27rz(M- 1)5^,01*^,0- LetM/N <e. 

k 1=0 

Then 

IIX2-X2II <7e 

Proof 

We have fjk^i = XkHk,i (look at (3), (4)), where Hk^i = Hk,i{A). Then = XkHk,i{-A). 
Now Hk,i{A) = i7fc/-A)exp(27rz(M - 1)A), \Sk,i - A| < 1/N, 2n{M - - A| < 7M/N 

and Lemma 2 follows. 



3.4 Complexity of a wizard 

It is readily seen that the above procedure for p = n requires applications of the initial 
transformation U, hence it cannot be used on purpose to predict its evolution. It turns that in 
general case the following theorem takes place. 

Theorem 1 It is impossible to simulate a wizard action with precision up to 0{1/N) using 
less than fl{N) conditional application of initial transformation Ucond- 
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Proof 

A lower bound for the time of quantum computation of the PARITY function is N/2 by the 
results ||1:''GGS|| and [ [BBCMWH - These constructions can be extended to the case when we can 
use more general oracle Ucond instead of U when evaluating PARITY. Assuming that a wizard 
action can be simulated using less than VL{N) application of Ucond, in view of the previous 
subsection we would be able to compute the function PARITY on a quantum computer in time 
less than N/2 (adding appropriate constant to the number n in order to enhance accuracy) 
which leads to the contradiction. Theorem 1 is proved. 

Call a unitary transformation U classical if it maps basic states to basic states. One can 
ask is there any way to predict a behavior of classical system on a quantum computer? This 
statement is independent of a form of spectrum. It turns that in general case a behavior of 
the bulk of classical systems in a short time frame is impossible on quantum computer even on 
one step ahead (look at |P^| ). Here short time means approximately 0{N^^'^)^ where N - the 
number of all states. 

But for one class of systems the quantum prediction is possible. This is the class of systems 
with sparse spectrum. 



4 Prediction of the evolution of the systems with sparse 
spectrum 

4.1 Case of precise eigenvalues 

What would happen if we use less number of qubits p < n instead of n in the ancilla? Simulation 
of the wizard will be shorter, but it will be simulated with the corresponding loss in precision 
and we obtain no prediction. 

Nevertheless, if the spectrum is sparse then prediction is possible. Indeed, suppose that we 
have classical algorithm h enhancing the accuracy of eigenvalues approximation. 

Namely, let p < n, M = 2^ and let h : {0, 1}^ — > {0, 1}" be integer function mapping 
rough frequency approximation up to 1/M to the more precise approximation up to 

Given the initial state ^ and time instant t we shall now describe the procedure of quick 
finding f/*^. The idea is simple: Repeat the procedure from above with only p ancillary qubits 
instead of n, and enhance accuracy by h. 

At first we assume for the simplicity that we have exact equalities uJk = ■ Here is the 
detailed description of the algorithm. 

First we apply Walsh-Hadamard transform to the ancilla and obtain 

N-1M~1 



k=0 s=0 

Then we apply U^gg. Since s < p this operation requires p conditional applications of U. 
This gives the state 

N-lM-l N-IM~1 

Xi = -Tfj E E XkU'\^k,s) = -j= E E exp(27rza;fcs)xfc|$fc,s). 
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Now application of QFTj^ to the second register yields: 

N-l 

X2= Yl ^k\4>k,'^k) 
k=0 

where O.Wk is approximation of eigenvalue Uk in within 1/M. 

This is the point when we will use that the spectrum is sparse. Add one more register with 

N-l 

n qubits initialized by zeroes: J2 a;A;|$fc, c^fc, 0"). Apply the unitary version of the algorithm 

fc=0 

h enhancing accuracy to two ancillary registers: \a,b) — > |a, 60/i(a)). This gives the state 

iV-l 

J2 Xk\^k, ^ki h{Cbk))- Now we turn every state of the form cD^, I) to the angle 27r -Q.l-t and 

k=0 

N-l 

obtain the state Xfccxp(27ri • OJ • t)\(^k,^k,l) ■ Here i is a given time instant. Now apply 

fe=0 

unitary version of h which cleans up n last ancillary qubits and then discard them: 

Af-l 

a;fcexp(27ru;fct)|$fc,'^fc), 

fc=0 

and again QFT^i^ to the p ancillary qubits: 

N-lM-l 

-j= Y Z] a;ikexp(27riu;fei)exp(-27ria;fes)|$jk,s), 

V M jk=o s=Q 

then again Ufg^ and Walsh- Hadamard transform to the ancilla. The result will be: 

N-l 

Xk cxp(27rzcjfet)|<l>fe, 0^). The wanted state is in the main register now. 

Note that all the work with enhancing an accuracy here is not mandatory because ujk = oj^ 
exactly. Why this scheme doesn't work for the case when Cj^ are only approximations of 
the true eigenvalues ujk ? The point is that X2 will not have so simple form and in addition 
multiplication frequencies by t will cause big error for t — 0{N). Thus for the general case 
more refined algorithm is needed. 



4.2 General case 

We have numbers n and p. Choose some q : p < q < n and put L = 2'^. Without loss in 
generality we can extend h to the mapping h : {0, 1}^ — > {0, 1}" so that if lj^ = {0.x)g then 

Our initial state has the form 

N-l 
k=0 

1. Apply QFT^U^ggWH to the initial state to get Xq,2- 
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2. Add one more register with n qubits initialized by zeroes: 

L-l 

yk,i\^k, 0^) and apply the unitary version of the algorithm h enhancing accuracy to two 

k 1=0 

L-l 

ancillary registers: Xs = E E yk,i\^k,l, h{l)) . 

k 1=0 

3. Turn every state of the form /, h{l)) to the angle 27r ■ 0.h{l) ■ t and obtain the state 

L-l 

X4 = E E yk,i cxp(27rz ■ 0.h{l) ■ t)|$fc, /, h{l)). Here t is given time instant. 

k 1=0 ' 

4. Put 61 1 = {0.h{l))n - (0./)g. Turn every state of the form to the angle 
— 27r(L — l)Sl i. Denote the result by X5- 

5. Now apply unitary version of h which cleans up n last ancillary qubits and then discard 
them: 

X6 = E E'l/Mexp(27ri • 0.h{l) ■ t) exp(-27ri(L - l)5li)\^k,l)- 

k 1=0 

6. Apply to xe the transformation WH U^g^ QFT^, then observe and discard ancilla. 

The time of this algorithm is 2L + 2iQFT^ + 2t/i + w where ^qft^ is the time of Fourier 
transform, th is the time of enhancing accuracy and w is the time of rotations. 



Lemma 3 Denote the result of algorithm 1-6 by Vq,t- For any 5 > for any p, q, n, t altering 
so that q — p + c, L — 2'^ , t : < t < CN, where c — log2(|); p > c, C — we shall have 



l|U,,t-U*|| <5 

Proof 

Introduce the simplifying notations: 
X = E E'l/Mexp(27riu;fci)|$fe,Z,/i(0>, X = QFT^^ {^ieg)~' WHexp(27ria;fei) E^^^l^fe, 0^), 

k 1=0 k 
L-l 

X' = E E yk,iexp{27riujkt) exp(-27ri(L - l)6li)\^k,l) ■ 

k 1=0 

By Lemma 1 ||x4 — xll < f • The passages Xi — ^ Xe si-nd x — ^ x' fulfilled by the same 
unitary transformation, which preserves lengths. Consequently, ||X6 ~ x'll < f • By Lemma 2 
llx — x'll < f- Then triangle inequality gives ||x6 ~ xll < ^- The passages Xe — ^ Ug^t and x — 
[/* are fulfilled by the same unitary transformation WH U^g^ QFT^. Hence \\U* — Uq^tW < 
Lemma 3 is proved. 

By Lemma 3 this algorithm gives the prediction of state in time instant 0{N) in time 
0{M) if classical algorithm enhancing accuracy obtain eigenvalues in negligible time. Thus if 
it is possible to enhance accuracy of eigenvalues from ei to e the speedup will be ^^-5^^, where 
p is a fidelity, < p < 1. The result can be formulated as 

Theorem 2 Given a Hamiltonian of system with sparse spectrum and algorithm enhancing the 
accuracy of eigenvalues from Si to e and a fidelity p, a state of the system in the time instant 
can he obtained in time ^i^p^^ ''^'ith this fidelity. 
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4.3 Generalizations 



4.3.1 Sparse areas of spectrum 

A spectrum of real system like molecule typically contains strips where spectrum is continuos 
separated by gaps where energy levels are absent at all. Let w be the maximal width of strips 
and g be the minimal width of gapes. 

Regard the times of evolution less than 1/w (here we assume the system of units where 
Plank constant is unit). For these times the width of strips is negligible and we can assume that 
El = 1/(7 in Theorem 2. So in this situation e = w and we obtain the following generalization 
of Theorem 2. 

Ift = o{l/w) we can predict the state U(t) of primary system in time tpre = 0{t^). 

Yet more generalization can be obtained if we consider the spectrum which is not sparse at 
all but has sparse areas. Here if the initial state of primary system is concentrated in sparse area, 
the prediction is possible with such probability which is equal to the degree of concentration. 



4.3.2 Travelling in a time: predictions and restoration of a history 

A wizard transformation may be used not only for predictions a future but for restoration of 
a history as well. If we replace the time t by —t in predicting procedure then we obtain the 
state of primary system in time instant —t which means the restoration of a history. Perform 
procedure from the section 4.2 with —t instead of t, where < t < CN. This makes possible 
to obtain the state of primary system in time of order —N in the past. Thus we conclude the 
following generalization of Theorem 2. 

Theorem 3 Given a Hamiltonian of some system, its states in time instant T = in future 
as well as in time instant —T in the past can be obtained in time with fidelity p provided 

each eigenvalue can be quickly calculated with precision up to e given its approximation with 
accuracy si. 

Note that for restoration of a history we don't need a sparse spectrum, may be Si = e. 
This is surprising because the algorithm uses only operator U when the natural way to obtain 
a state from the past is to use an inverse transformation and there is no evident way to 
simulate the action of the inverse transformation by means of U. 



4.4 Examples 

4.4.1 Shor factoring algorithm 

Formula (2) for a wizard simulation can be considered as a natural generalization for the 
quantum part of Shor factoring algorithm ( |ph|| ). Let a unitary operator U in (2) has the form 
U|x) — > \ax (mod g)), where ax is a multiplication of integers, (a, g) = 1. Then eigenvectors 

r-l 

of U have the form -j- eiqi{—2Txikj/r)\a^ (mod q)) and its eigenvalues are exp(27rij/r) where 

i=o 

r is a period of a (minimal integer such that = 1 (mod q)). If we apply an operator (2) 
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with this U to zero initial state and observe the second register then by Lemma 1 we obtain 
an approximation of a number j /r up to 0(1/ N) with high probabihty. Using this procedure 
sequentially we can restore the value r of a period. This procedure was used for factoring by 
Shor. 

Here we use a wizard simulation (2) as a technical element of the algorithm 4.2. Alter- 
natively, one can try to use Kitaev's method instead of QFT in (2). Namely, consider one 
controlling qubit and use Hadamard transform after conditional application of U and repeat 
this with numerous controlling qubits. Then by some auxiliary transformation it is possible 
to extract eigenvalues in ancillary register (look at ||Ki|| for details). But then we must get rid 
of revealed eigenvalues like in the point 6 of the algorithm and it faces a little difficulty when 
using Hadamard transform instead of QFT in (2), because we have not Uj^^^ and must always 
manage with only Ucond- 

4.4.2 Grover search algorithm 

Consider Grover algorithm for the fast quantum search: U{t)\0) = (/Q/a)*|0), where denotes 
inversion of the state b: Ii,\b) = —\b), Ib\s) = \s) for {b\s) = 0. It was shown in | |Gr2| ] that 
if |(0|a)| = 0(1/ VN) then U{h) ^ \a) for ti = 0{VN) independently of 0, where N is the 
dimension of main space. What will happen if we apply Theorem 2 for U as a primary evolution? 
The problem in fact will be two dimensional and the minimal gap between eigenvalues will be 
of order 1/\/N. Hence, by Theorem 2 we can predict the states U{t) for t = 0{N) in time of 
order y/N. As for t = 0{\fN) we obtain no additional speedup. 

4.5 About practical implementations 

The two main components of the algorithm are conditional iterations of a primary transforma- 
tion U and QFT. 

The main difficulty for the practical implementation of this method is in conditional itera- 
tions of U. Given only a primary device realizing U one cannot immediately fulfil this iterations 
because it requires quantum control on the number of iterations. The solution may be following. 
Decomposition of a primary device into elementary parts that can be included to the quantumly 
controlled circuit and realize the conditional iterations of U for all parts simultaneously. All 
controlling qubits should be used in entangled state as a Shroedinger cat. Given a quantum 
gate array computing U we can easily construct a new gate array computing X^cond- To put it 
otherwise a conditional application of U is possible through the control on microscopic level 
(this bears a resemblance with the control in living cells). 

This scheme can be reformulated by means of analogous quantum computing if we consider 
QFT on an ancillary register as a passage to the canonically conjugate magnitude. Say, if we 
use a value of coordinates in operations with a register then canonically conjugate will be the 
corresponding impulse. The passage from the coordinate representation of a wave function to 
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the impulse representation in one- dimensional case can be defined as 

+00 

4'ip) = J ex.p{—ipx / h)ip{x)dx 
—00 

where a probability to obtain impulse in a segment {p,p + dp) is \(j){p)\'^dp/27ih. Assume that 
we have one particle which can be located in M points of the form x = 0, j^, jj, . . . , 
Then the coordinate quantum space for one particle will be M dimensional. Consider the 
corresponding integral sum for (f){p) in the system of units where Plank constant h is one. This 
integral sum will be just the sum in the definition of QFT where x plays a role oi s/M and 
p/2'K plays a role of /. 

Then the main algorithm acquires the following general form. 

a) . Primary evolution quantumly controlled by a magnitude containing in the properly 
prepared ancillary register. 

b) . Simple actions depending on the canonically conjugate magnitude. 

c) . Repetition of a). 

Procedure of such a kind can predict a behavior of a primary system with sparse spectrum 
and restore the state of arbitrary system in the past. This scheme seems to be very simple 
and would be interesting to find its natural physical analog. The good starting point here may 
be the comparison between spectral features of the known systems and their functions and 
complexity. 

5 Conclusion 

Formulate again the main result: the method is presented which makes possible to obtain states 
of a primary quantum system earlier than they appear naturally in its evolution provided the 
spectrum of system is sparse. This speedup will be the more if the gaps between continuous 
strips of spectrum increase comparatively to the width of strips. This method can be applied 
also for restoration of states of the arbitrary primary system in the past. This method yields a 
speedup and it can make possible to fit into the time when coherent states exist and thus fight 
decoherence in quantum computations. 

6 Acknowledgements 

I am grateful to Lov Grover for his kind invitation to Bell Labs, discussions and attention to 
my work. I also thank Peter Hoyer and David DiVincenzo for fast replies and useful comments. 

References 

[AL] D.S.Abrams, S.Lloyd, A quantum algorithm providing exponential speed increase for 
finding eigenvalues and eigenvectors lanl e-print \quant-ph/ 980707^ 



11 



[BBCMW] R. Beals, H. Buhrman, R.Cleve, M. Mosca, R. De Wolf, Tight Quantum Bounds 
by Polynomials, lanl e-print \quant-ph/980204^ 

[FGGS] E. Farlii, J. Goldstone, S. Gutmann, M. Sipser, A limit on the Speed of Quantum 
Computation in Determining Parity, Phys. Rev. Lett., 81, 1998, 5442-5444, lanl e-print 
\quant-ph/98U2U45i 

[Gr] L. K. Grover, A fast quantum mechanical algorithm for database search. Proceedings, 
STOC 1996, 212-219. Philadelphia PA USA, lanl e-print ^uant-ph/96U5U4B[ 

[Gr2] L.K. Grover, Rapid sampling through quantum computing, lanl e-print \guant-ph/99 12001 

[Ki] A. Yu. Kitaev, Quantum measurements and the Abelian Stabilizer Problem lanl e-print 
\quant-ph/951102(\ 

[Oz] Y.I.Ozhigov, Quantum computers speed up classical with probability zero Chaos, Solitons 
and Fractals, 10(10), 1999, 1707-1714, lanl e-print \qumt-ph/ 980306^ 

[Sh] P. W. Shor, Polynomial-time algorithms for prime factorization and discrete logarithms 



on a quantum computer, SIAM J. Comp. 1997, 26, No. 5, 1484-1509, lanl e-print \quant- 
ph/950802^ 



12 



